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Abstract. The steady problem resulting from a mixture of two distinct fluids 
of power-law type is analyzed in this work. Mathematically, the problem 
results from the superposition of two power laws, one for a constant power- 
law index with other for a variable one. For the associated boundary-value 
problem, we prove the existence of very weak solutions, provided the variable 
power-law index is bounded from above by the constant one. This result 
requires the lowest possible assumptions on the variable power-law index and, 
as a particular case, extends the existence result by Ladyzhenskaya 1171 to the 
case of a variable exponent and for all zones of the pseudoplastic region. In a 
distinct result, we extend a classical theorem on the existence of weak solutions 
to the case of our problem. 
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1. Introduction 

In Fluid Mechanics the constitutive equation relates the stress on a fluid element 
to other fluid quantities by the relation 

T = -pi + S , 

where T is the Cauchy stress tensor, p stands for the pressure, I is the unit tensor, 
^ is the dynamic viscosity and S is the deviatoric part of the stress tensor. The 
simplest model of Fluid Mechanics is the Newtonian fluid which by definition is a 
fluid that continues to flow, regardless of the forces acting on it. For these fluids, 
the constitutive equation is the Stokes law and therefore the deviatoric part of the 
stress tensor is given by 

S = 2//D, D = D(u)^i(Vu + Vu^) , 

where D is the rate of strain tensor and u is the velocity field. Although the New- 
tonian fluid model captures the characteristics of many fluids as water solutions, 
gasoline, vegetal and mineral oils, this model is quite inadequate for describing the 
complex rheological behavior of many other fluids. In this case are all the fluids 
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in which the apparent fluid viscosity /i decreases or increases with the shear rate 
history |D|. These fluids form a large class called non- Newtonian, or generalized 
Newtonian fluids, and are divided into pseudoplastic and dilatant fluids. In pseudo- 
plastic fluids the viscosity gradually decreases with increasing shear rate and, due 
to this, they are often called shear thinning fluids. Shear thickening fluids is the 
other name found in the literature for dilatant fluids and are fluids in which the 
viscosity increases with the shear rate. Examples of pseudoplastic fluids are milk 
fluids, varnishes, shampoo and blood fluids, whereas polar ice, glaciers, volcano 
lava and sand are all examples of dilatant fluids. There are also some another class, 
called Bingham fluids, that are similar to pseudoplastic fluids, but they exhibit a 
yield point. Examples of this fluids are drilling muds used in petroleum industry, 
toothpaste and face creams. The existence of a yield point means fluid flow is pre- 
vented below a critical stress level, but flow occurs when the critical stress level is 
exceeded. During the 20th century a large number of models have been proposed 
in the literature to model all types of non-Newtonian fluids under diverse flow con- 
ditions. However, it was only since the pioneer work done by Oldroyd |20| . during 
the 1950's, that were established some guiding principles to provide a constructive 
framework for the formulation of constitutive equations to the Cauchy stress tensor 
(see Barnes et al. [S]). The power-law, or Ostwald-de Waele model, is one of the 
simplest non-Newtonian fluid models then introduced and can be defined by the 
relation 



where 7 > 1 characterizes the flow behavior and is usually called the power-law in- 
dex. The power-law is often used to model pseudoplastic fluids though it can also be 
used for modeling dilatant fluids or for Newtonian fluids (see e.g. Schowalter |27|): 



There are also some fluids that cannot be cataloged into a single class of non- 
Newtonian fluids. These fluids can go, for instance, from the consistency of a liquid 
to that of a gel, and back, with response times on the order of milliseconds. In this 
case are the electrorheological fluids whose rheological properties are controllable 
through the application of an electric field, showing useful and special function 
with the effect of reversibility. The best example of electrorheological fluids are 
suspensions dispersed with some polymeric colloids which show trembling shear 
behavior under an electric field (see e.g. Hao [H]). Due to this, in the sequel, by a 
trembling fluid we mean a fluid with a variable power-law index, whereas sustaining 
fluid is the name we shall use for a fluid with a constant power-law index. In this 
way, trembling fluids can also be medelled by the power-law but with the 

significant difference that now the power-law index 7 may vary with other quantity 
under study or, in the simplest cases, with the space and time position (see e.g. 
Rajagopal and Riizicka [24]). 

In this article we will study the mathematical problem for a steady motion of 
a generalized fluid contained in a bounded domain C R^, N > 2, with the 
boundary denoted by dfl. We assume the motion is described by the following 



(1.1) 



S = /i|D|T-2D, 



power- 




if 7 = 1 

if 1< 7 < 2 



if 7 = 2 
if 7 > 2 . 
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boundary-value problem for the generalized Navier-Stokes equations: 
(1.2) divu = in fJ; 



(1.3) 



div(u ® u) = f 



Vp + divS in fi; 



(1.4) 



u 







on dfl; 



where, p stands, now, for the pressure divided by the constant density and f is the 
external forces field. We assume the dependence of the deviatoric stress tensor S on 
the space variable x and on the strain rate tensor D is given by the superposition 
of two different power laws, one for a sustaining fluid with other for a trembling 
one: 



Here, and /i2 are positive constants related with the fluid viscosity, and 7 and 
q are the sustaining and trembling power-law indexes, respectively: 7 is a constant 
and q depends on the space variable. The specific physical problem we are interested 
in and where (|f .5p can be potentially useful is the case of a dilute suspension of 
an electrorheological material in a dilatant fluid. The object of superposition of 
generalized fluids is to produce flow patterns similar to those of practical interest. 
Many systems, among them polymer solutions and emulsions, behave, in the dilute 
regime, as superposed fluids in the above sense. The best example are polymer 
solutions in which the polymer segments tend to repel each other, since they prefer 
contact the solvent molecules rather then among themselves (see e.g. Oswald \'22\). 
Moreover, superposition of fluids is justified, in the light of theoretical mechanics, 
as a powerful tool to replace the Boltzman superposition principle in the case of 
materials with nonlinear behavior (see e.g. Dealy (TOl). The constitutive relation 
(|1.5p can be used to model many other generalized fluids as follows: 



In particular, by making 7 = 2 in ()f .5p . we obtain a generalization of the Sisko 
model to the trembling fluids. The Sisko model has been checked experimentally to 
fit accurately the viscosity data of various commercial greases made from petroleum 
oils with one of the standard thickening agents such as calcium fatty acid, lithium 
hydroxy stearate, sodium tallow or hydrophobic silica over a wide range of shear 
rate (see Sisko l29j). Letting also 7 = 2, we recover a specific Carreau-Yassuda 
model which is very often used to describe blood flows (see e.g. Carreau et al. [8]): 



In the specific example we want to address, //q stands for the zero-shear viscosity, 
the infinite-shear viscosity yb^o is zero, where A is a relaxation time, the power-law 
index is absent, i.e. n = 2, and a stands for a variable shape parameter. 

The outline of this work is the following. In Section [1] not only we presented 
the problem that we shall study in this work, but also we have given a physical 
motivation for doing so. The main notation used throughout the text and some 
auxiliary results are presented in Section [5] Section [3] is devoted to review the main 
existence results for some particular cases of the problem (|1.2p - (|f .5p . In Section[3]we 
define the notion of solutions we shall consider and we state two different existence 



(1.5) 




sustaining power-law if /^i > and /i2 = 
trembling power-law if /ii = and /i2 > . 



S = A'oo + (mo - A^oo) (1 + A|Dr) | — 
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results: Theorems 14.11 and 14.21 From Section [5] to Section [10] we shall prove the 
main result of this work: Theorem |4JJ 



2. Preliminaries 

The notation used in this work is largely standard in Mathematical Fluid Me- 
chanics (see e.g. Lions [TH]). In this article, the notations il or lu stand always for 
a domain, i.e., a connected open subset of E^, TV > 1. Given fc G IN, we denote 
by C^{il) the space of all fc-differentiable functions in il. By Cg°(n) we denote 
the space of all infinity-differentiable functions with compact support in ft. In the 
context of distributions, the space Cjf (f2) is denoted by 'D{il) instead. The space 
of distributions over 2?(r2) is denoted by 2?'(i7). If X is a generic Banach space, 
its dual space is denoted by X'. Let I < q < oo and ft C R^, with > 1, be 
a domain. We use the classical Lebesgue spaces L^(17), whose norm is denoted by 
II ■ I|l<!(si)- For any nonnegative fc, W'^''^(fi) denotes the Sobolev space of all functions 
u £ L'^{il) such that the weak derivatives D"u exist, in the generalized sense, and 
are in L*(51) for any multi-index a such that < |a| < k. In particular, W^'°°{il) 
stands for the space of Lipschitz functions. The norm in W''^''{n) is denoted by 
II • llw^.^fJ)- We define Wo'''(r2) as the closure of C^ifl) in W'=^'?(rj). For the dual 
space of Wg'^(r2), we use the identity (WQ'^(fi))' = W"*^'^ (il), up to an isometric 
isomorphism. We shall distinguish (second-order) tensor- valued and vector- valued 
space functions from scalar-valued ones by using boldface letters. Although we are 
going to use the same notation for tensor-valued and vector-valued space functions, 
the distinction between these spaces will always be clear from the exposition. 

We denote by V{fl) the set of all measurable functions g : J7 — >■ [1, c»] and define 

:= cssinf g(a;), q^ := esssupg(a;). 

Given q G V{n), we denote by L«( )(rj) the space of all measurable functions / in 
n such that its semimodular is finite: 



(2.1) ^,(.)(/):= / |/(x)r(-'dx< 



oo. 



The space L«(-) (fi) is called Orhcz space and is also known by Lebesgue space with 
variable exponent. Equipped with the norm 

(2.2) II/IIl.()(o) :=inf|AJ>0:A,(.) (^0 <1 

L'(')(ri) becomes a Banach space. Note that the infimum in (|2.2p is attained if 
^q{-)if) > 0- If 1^ < L'i'( '(ri) is separable and the space Cjf (fi) is dense in 
L'J(-)(f]). Moreover, if 

(2.3) 1< q^ < q+ < oo, 

L'^^'^n) is reflexive. One problem in Orlicz spaces, is the relation between the 



semimodular (|2.1|) and the norm (|2.2|) . If (|2.3p is satisfied, one can shows that 

(2-4) ii/ii^;,(,,)-i<^(.)(./)<ii/ii^:,,(n) + i- 

In Orlicz spaces, there holds a version of Holder's inequality, called generalized 
Holder's inequality. 
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Given q <E 'P(f^), the Orlicz-Sobolev space is defined as: 

:= |/ e Li^-\n) : D"/ e L''( )(f]), < |a| < l} . 

In the hterature, this space is also known by Sobolev space with variable expo- 
nent. In VF^'''^'^ (rj) is defined a semimodular and the correspondent induced norm 
analogously as in (|2.ip - (|2.2p . which are equivalent, respectively, to 

^.,(.)(/) :=A,(.)(/) + ^,(.)(|V/|) 

and 

■= II/IIli( )(J2) + l|V/||L9( )(n) • 

For this norm, is a Banach space, which becomes separable and reflexive 

in the same conditions as L''(')(ri). The Orlicz-Sobolev space with zero boundary 
values is defined by: 

:= {/ g w^M-){n) : supp / cc n} "■""-^""e^) . 

In contrast to the case of classical Sobolev spaces, the set Cj^(r2) is not necessarily 
dense in W^'^^'^n) - the closure of C^{n) in W^-'^'-'^n) is strictly contained in 
Wq'''' ■'(SI). The equality holds only if q is globally log-Holder continuous, i.e., if 
exist positive constants C\, C2 and qoo such that 

(2.5) |g(x)-g(y)| < f ; |g(x)-g^|< i n ^ x, y e f]. 

ln(e + l/|x - y|) ln(e + |x|) 

For a thorough analysis on Orlicz and Orlicz-Sobolev spaces, we address the reader 
to the monograph by Diening et al. |11| . 



3. Historical background 

To the best of our knowledge, problem (|1.2|) - (|1.5p is new, but there is an extensive 
literature on some particular cases of it. Here we shall be concerned only with the 
existence results for the problem (|1.2p - (|1.5p or its simplifications. 

In the case of /xi ^ and 112 = 0, the main achievements on the existence results 
for this problem were done by Lions [19], Frehse et al. |131fT4] and Riizicka [25]. But 
even earlier than the first of these authors, the mathematical analysis of the prob- 
lem (|1.2p - (|1.5p with /ii 7^ 0, but q = 2 and A'' = 3, was performed by Ladyzhenskaya 
in the work [T7]. It should be noted that Ladyzhenskaya's problem [T7] is, in fact, 
the problem considered by Sisko [29], more or less 10 years before, in rheological 
experiments. Surprisingly is that, in the literature, there is no evidence that La- 
dyzhenskaya knew Sisko's work. Moreover, while Sisko was working on experiments 
with various commercial grease flows, Ladyzhenskaya studied this problem moti- 
vated by the important issue of uniqueness for the classical Navier-Stokes problem. 
Problem (|1.2p - (|1.5p with /ii 7^ and fi2 ^ was consider by Lions [19] using the 
same approach of Ladyzhenskaya ^7\. In both works [T7] and [19], and under the 
assumption that f £ V' and 



(3.1) 



3iV 

- N + 2' 
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the authors have proved the existence of weak solutions in the class V-y satisfying 
to the energy relation 

(3.2) / (S(D(u))-u«)u) :D((/3)(ix= / i-ipdx 
Jn Jn 

for all ip G V-y. Here the notation V'^ stands for the dual space of V-y and 

(3.3) V ■.= {veC^{n)■.div^r = 0}, 

(3.4) := closure of V in W^'''{n). 

The proofs in |17[ I19j use the theory of monotone operators together with compact- 
ness arguments, and the lower bound p.ip results from controlling the boundedness 
of the convective term u (g) u : D((/3) in L^(J7) for u and (p belonging to V^. More 
or less 30 years later the results of [171 IT^ were improved in [131 [2S] for values of 7 

such that 



N + 1 

Assuming that f G L''' (fi) in [T3] and f e V^^ in ^23), the authors have established 
existence results in the same class of [ITl [19] satisfying to the energy relation 
p.2p . but for all (y3 G V. In order to distinguish the weak solutions found in \Tf\ [T9] 
from these established in |13[ I25| , we shall denote the later as very weak solutions. 
An important feature of the works [13) [25] , is that there the space of test functions 
allows us to consider test functions with more regularity. Due to that, to control 
the boundedness of div(u (E) u) ■ ip in L^(ri) for u in and ip in V, led us to 
a lower bound for 7 (/. ()3.5I1 ). Although this important difference, the proofs of 
[T3] and |25] also use the theory of monotone operators together with compactness 
arguments and they differ only on a slight different application of the L°°-truncation 
method. A little bit later, in [T3] the authors improved theirs own result [T3] for 

2N 

(3-6) 7 > tF^- 

^ ' ' N + 2 

The new feature was the application of the Lipschitz-truncation method. By this 
method, the authors could use all the regularity they needed for the test function 
and to worry only about the boundedness control of u (g) u in (il) for u in 
and p in V, which holds for 7 > . The strict inequality in (j3.6[) results from 
the validity of the compact imbedding V-^ L^(il), which is of fundamental 

use, not only in [14) . but also in all aforementioned works. 

The analysis of the problem (|1.2p - (|1.5p with fii ~ Q and /i2 7^ started, from 
the Mathematical Fluid Mechanics viewpoint, with the works on electrorheological 
fluids by Riizicka [25] and on thermorheological fluids by Antontsev et al. [2l [3l [4] . 
The main existence results for the problem (|1.2l) - (|1.5p . with fii = and /i2 ^ 0, are 
due to Ruzicka j^S], Huber |TB] and Diening et al. [T^]. The first existence result 
for the problem (|1.2p - p.5p . with yUi = and /i2 7^ 0, is an immediate consequence of 
the same result to the corresponding electrorheological problem established in [2B] . 
In fact, proceeding as in the proof of [26] Theorem 3.2.4], one can easily proves the 
existence of weak solutions to our problem in the following class 

(3.7) '^q{-) closure of V in the ||D(v)||L<j(.)(f^)- norm, 
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and satisfying to the energy relation p.2p for all f E Wq(.). Moreover, the result 
holds for q G V{fl) by assuming that 

(3.8) 1 < a := cssinf ^(x) < g(x) < esssup(7(x) := /3 < oo, 

f € and a satisfies to (|3.ip . The proof here follows the same approach of |171 [T^ 
and uses the fact that Wq(.) is imbedded into Vq. In [TC] the existence result [55] 
was improved for the case of a satisfying to p.5p . Here, under the assumptions 
that f £ (Wq''^' -' (n))' and (|2.5p holds, is proved the existence of weak solutions in 
the class 

(3.9) V^f.) := closure of V in Wi-9( ) (fl) 

satisfying to the energy relation (|3.2[) for all (/3 G V. The proof there combines 
a generalization of Bogowskii [6] results on divergence problems in Orlicz-Sobolev 
spaces with the approach followed in |131 [25] for the case of a constant q. Finally, 
under the same assumptions of (TH], it is proved in [T2] a more complete existence 
result in the class V^^.-). This result holds for a satisfying to p.6p and the solutions 
satisfies to the energy relation 

(3.10) / (S(D(u)) - u) : D((^)(ix = / pdiv^(ix+ / i-ipdx 
Jn Jn Jn 

for all ip £ WQ'°°(f2). The proof follows the same approach of the result for constant 
q by the same authors and uses results on Lipschitz truncations of functions in 
Orlicz-Sobolev spaces performed by the authors still in |12j . 

A crucial assumption in the works |12[ I16| is that the variable exponent q must 
be globally log-Holder continuous. Recall that q is globally log-Holder continuous, 
if q is locally log-Holder continuous (/. p. 5^ 1 ) and if there exist constants C2 and 
qoo such that ()2.5p 9 holds. This assumption gets into the proof in order to use the 
denseness of Cjf^(f2) in W^'^' -'(f2) and, due to that, it is possible to look for the 
solutions in the class V^^.-j instead of W^^.). However, though for many physical 
problems of trembling fluids, the variable power-law index q satisfies the log-Holder 
continuity property (j2.5p . there are some mathematical studies that suggest the 
non sufficiency of the condition (|2.5p . For a discussion on this important issue, see 
Diening et al. Section 5.1] and the references cited therein. In the current work 
we shall establish an existence result for the complete problem (|1.2p - (|1.5p without 
requiring the variable exponent q is globally log-Holder continuous. 

4. Weak formulation 

In order to introduce the notion of weak solutions we shall consider in this 
work, let us recall the well-known function spaces of Mathematical Fluid Mechanics 
defined at (|3.3p - p.4p . Due to the presence of the variable exponent q{-) in the 
structure of the deviatoric tensor S, we need to consider the weak solutions to the 
problem (|1.2p - p.5p in some Orlicz-Sobolev space. As we already pointed out at 
Section [2l an important limitation of Orlicz-Sobolev spaces is that, without any 
extra condition on the variable exponent q, but q G Vi^l) satisfying to p.Sp . the 
set Cg°(J7) is not necessarily dense in Wq'''^ \^)- For this reason, we shall look 
for our solutions in the function space Wg(.) defined in p.7p instead of the one 
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considered in p.9p . It is a easy task to verify the space Wq(.) satisfies to tlie 
following imbeddings: 

(4.1) ^ W,(.) ^ V„ . 

Moreover. Wg(.) is a closed subspace of Vq and therefore it is a reflexive and 
separable Banach space for the norm 

l|v||w,(., := ||D(v)||L,(-)(n). 

Definition 4.1. Let fl be a bounded domain of M,^ , with N > 2. Assume that 
f G Li^{fl), J is a constant such that 1 < 7 < oo and q G 'P{^) is a variable 
exponent satisfying to iS.S]) . A vector field u is a (very) weak solution to the 
problem {UP-fTg), if: 

(1) ue W,(.)nv^; 

(2) For every ip G Wg(.) n (For every (p Cz V) 

J (A'i|D(u)r"' + M2|D(u)|«(")-2-u0u) :D(vj)dx = ^ f -(^dx. 

Remark 4.1. Note that ifj>l3, then V-y ^ '^q(-) '^'^'^ therefore it is enough to 
look for weak solutions in the class Wq(.). 

As we had mentioned by the end of last section, in this work we are mainly 
focused on existence results for the problem (|1.2p - (|1.5p . or for some of its simpli- 
fications, without invoke the log-Holder continuity property (j2.5p on the variable 
exponent q. We this in mind, an existence result for the problem that can be 
adapted from already known results is written in the following theorem. 

Theorem 4.1. Let il be a bounded domain in , N > 2, with a Lipschitz- 
continuous boundary dfl. Assume that 1 < 7 < 00, q e 'P(r2) satisfies to US. 8\} and 
f e (V^nW^(.))'. Then, if 

mm {7, a} > -j^—^ , 
there exists a weak solution to the problem 1773)) in the sense of Definition\4-.l\ 



Proof. The proof combines the result addressed in |19l Remark 2.5.5] for the con- 
stant power-law index 7 with the existence result of |26[ Theorem 3.2.4] for the 
variable power-law index q. We observe that, since (|4.ip holds, \26\ Theorem 3.2.4] 
is a simple extension of the result established in [19]. Therefore, we can say that 
the proof Theorem 14.11 follows easily from the result invoked in [19] . □ 



Remark 4.2. In the particular case 0/7 = 2, Theorem \4.1\ extends the existence 
result established by Ladyzhenskaya in |17j to the case of a variable exponent q, 
without invoking the log-Holder continuity property i2.5\) . Moreover, since V^j.j C 
Wg(.), this result is obtained in a larger class than the ones which follow from the 
works [T2l[T6l. 



The main result of this work is written in what follows. We establish here the 
existence of weak solutions for the problem p.2p - p.5p with a variable q depending 
on X and without any other restriction on q but (|3.8p . Our existence result will be 
valid for the constant exponent 7 and the variable one related via the relation (|4.2p 
below. 
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Theorem 4.2. Let fl be a bounded domain in M,^ , N > 2. Assume that q E 'P{^1) 
satisfies to i3.8]) and f E (V^ n Wg(.)) . Then, if for any 5 > 

(4.2) ^>niax|^^ + <5,/3| , 

there exists a very weak solution to the problem i [j.i?|) -i f7T3|) in the sense of Defini- 
tion \4-l\ 

The proof of Theorem 14 . 21 wiU be the aim of the next sections. In order to make 
the proof as transparent as possible, we shall assume that 

(4.3) f=-divF, FeM,^,„, FeL«'(-)(r!), 

where is the vector space of all symmetric N x N matrices, which is equipped 

with the scalar product A : B and norm |A| = VA : A. Note that assumption 
(|4.3p does not affect at all the extent of Theorem 14.21 because f = — div F and 
F e L«'(-)(f2) implies that f e W^^.^ and 7 > /3 implies W^^.^ ^ (V^ n W,(.))'. 
The assumption F € M^j^^ is made in order to avoid unnecessary calculus. But, 
before we get into the details of the proof, let us make a few comments. Firstly, 
we note that, contrary to Theorem 14.11 here we do not assume any regularity on 
the boundary dCl, besides the one following from the boundedness of fl. Then, we 
observe that, as in Theorem l4.1[ the only assumption on q, besides p.Sp and (|4.2p . is 
that q is a measurable function from fl into [1, 00]. Moreover, Remark l4.2l extends to 
this case of very weak solutions, now with 7 = 2 and q related via (|4.2p . This result 
is not mentioned in the literature, because, when q is constant, Ladyzhenskaya's 
problem is, from the mathematical viewpoint, equivalent to the problem treated 
in Lions |19| . This equivalence still remains in the case of a variable exponent q 
satisfying to the log-Holder continuity property (|2.5p (/. Diening et al. [12j). But 
if this hypothesis is not required, this result is, to the best of our knowledge, new. 
We observe that also the existence result established in Theorem 14.21 is valid for a 
larger class of functions than the results provided by the works [T2t [T6] . 



5. The regularized problem 



Let $ G C°°([0,(X))) be a non-increasing function such that < <& < 1 in [0,oo), 
$ = 1 in [0, 1], 3> = in [2, 00) and < < 2. For e > 0, we set 

(5.1) $,(s) := $(es), se[0,cx)). 

We consider the following regularized problem: 



(5.2) 

(5.3) 
div(Uf 



(5.4) 



div Ue = in f2. 



u,$,(|u,|)) Vp, 



div 



/ii|D(u,)r-^ + ^2|D(u,)r(")-') D(u, 







on dVl. 



Note that we have introduced a regularization that allows us to control the convec- 
tive term. A vector function u,: € is a weak solution to the problem (|5.2p - (|5.4p . 
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if 



f [(//i|D(u,)r-' + M2|D(u,)r(")-')D(u,)-u,®u,$,(|u,|)] :D(^)dx 

= [ F : D{^)dx 
Jn 

for all e V. Note that p.8p implies ^ ) f is given in the form (|4.3p . 



(5.5) 



Proposition 5.1. Lef f/ie assumptions of Theorem^J^ be fulfilled. Then, for each 
e > 0, there exists a weak solution G to the problem I15.^) - (S^^ . In addition, 
every weak solution satisfies to the following energy equality: 

(5.6) / (Aii|D(ue)r + M2|D(u,)|«(''))dx= / F:D(u,)dx. 

Jn ^ ' Jn 

Proof. Let us set 

M := {n7 e H : < 1} 

and lets us consider the following system: 

(5.7) divw = in f], 
(5.8) 



-div 



Mi|D(u)r-2 ^ ^2|D(u)rW-2^ D(u)' 



=f - Vp 
- div(ra7 (g) w<^e{\'^\)) 

(5.9) v = Q on dVL. 

Observing that, due to (|4.2p . the continuous imbedding ^^{ri) ^ L'^( '(il) holds, 
we can use the theory of monotone operators (see e.g. Lions [191 Section 2.2]) to 
prove that for each vj € M , there exists a unique weak solution v S to the 
system ^1^-^^. 

As a consequence of the previous step, we can define a mapping 

(5.10) K : M ^ H 

such that to each w £ M associates a unique ~ ^v, where v G V-^ is the unique 
weak solution to the system (|5.7p - (|5.9p and ^ is a positive constant which will be 
defined later on. Testing formally (|5.8p by the unique weak solution v such that 
S,v := K(i37), with zu G A/, integrating over Q, using the inequalities of Young and 
Korn. and at last the definition of $£(|i37|). we achieve to 

Ci f (|D(^)r + |D(^)rW)dx<7i+72||HlL(n), 

(5.11) 

71 := C2 / iFp'dx, 72 := C3. 
Jn 

Then, setting ^ (C^/Cii^i + 72))"^, we can prove, using the inequalities of 
Sobolev and Korn, the generalized Holder's inequality, and then (|5.1ip . that 

(5.12) ||K(n7)||L.(o) <eC4 (^^|D(w)rWdx+l^ <CC4/Ci(7i+72) = 1 

for all w £ AI. This proves that K maps Af into itself. 

In order to prove the compactness of K, we observe that from (|5.1ip it follows 

(5.13) ||K(tn)|i^^ <C(7i+72) 
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for all ZD e M. Owing to the embedding \fl '■'''(^) '"^ L''' (Jl) and to the assumption 
(|4.3p . the right hand side of (|5.13|) is finite. Then, due to the compact imbedding 



L^(f^), valid for any 7 > K(Af) is relatively compact in L^(ri). 

To prove the continuity of K, we consider a sequence in M such that 

vjra w in L^(il) as m — > 00. 

By the relative compactness of K(Af) in L^(ri), there exists a subsequence Wmt, 
such that 

(5.14) K(ro„jj^) — > <; in L^(ri), as k ^ 00. 

From the definition of K, the functions defined by ^ = Tnii. K(n7m^), 
satisfy to 



(5.15) 



(/ii|D(z;™jr-2 + ^2|D(z;„jrW-2) D(^;„J : D(^)dx = 

Jn 



for all (/? e V. Passing to the limit in (|5.15p and using the Minty trick, we can 
prove that <^ = and therefore = K.{zu). Then, from (|5.14p . we conclude that 
K(-ci7m) — > K(n7) in L^(ri) as m 00, which proves the continuity of K. 

Now, applying Schauder's fixed point theorem, there exists a function e M 
such that K(u^) = and which is a weak solution to the problem (|5.2p - (|5.4p in fl. 

Finally, the energy relation (|5.6p follows by testing (|5.3p by a weak solution and 
integrating over f2, and observing that now the convective term is zero. □ 

6. Existence of approximative solutions 

Let Ug € be a weak solution to the problem (|5.2p - (|5.4p . From Proposition [5TT] 
(/. (|5.6p '). using the generalized Holder's inequality and the relation between the 
norm ||D(u£)||L,(.)(f2-) and the semimodular (D(u<;)) (/. (|2.4p ') together with 
Young's inequality, we can prove that 



(6.1) 



^(|D(u,)r + |D(u,)rW)dx<C, 



where, by the assumption (j4.3p . C is a positive constant and, very important, does 
not depend on e. Appealing again to the relation between the norm ||D(ue)||L<j( )(f2) 
and the semimodular A^(-.-)(D(ue)) and also to generalized Holder's inequality, we 
can prove from (|6.ip that 

(6.2) ||u,||v, < C\ 

(6.3) ||D(u,)||l,(.,(o) < C. 
Proceeding as for (|6.2p - (|6.3p . we can also prove that 

(6.4) |||D(u,)r-2D(u,)||Ly(o) <C, 

(6.5) |||D(u,)r('')-2D(u,)|lL,v,(o)<t^- 
Moreover, using (|4.2p . it also follows that 

(6.6) |l|D(u,)rW-2D(u,)|Ly(o)<C. 
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On the other hand, by using (|6.2|) and Sobolev's inequahty, we have 

(6.7) \\neh.'in)<C, 

where 7* denotes the Sobolev conjugate of 7. As a consequence of ()6.7p and due 
to the definition of $e (/. (|5J|l ). 

(6.8) ||u,®u,cl>,(|u,|)|l^^^^^ <C. 

Note that the constants in (|6.2l) - (|6.8p are distinct and do not depend on e. From 
(|6.2p . (|6.4p . (|6.6p and (|6.8p . there exists a sequence of positive numbers e,„ such 
that e„i — > 0, as m — > 00, and 

(6.9) Ue„^ — > u weakly in V-^, as m ^ 00, 

(6.10) \B{u,Jp~'^'D{u,^) ^ Si weakly in L'^'(17), as m^oo. 

(6.11) |D(u,„J|''W-2D(ue,J ^ S2 weakly in L'''(17), as m 00. 

(6.12) (g) Ue^$g^(|ue,^|) ^ G weakly in Lt" (fi), as to od. 



(6.16) 



Now we observe that, due to (|6.9p . the application of Sobolev's compact imbed- 
ding theorem implies 

(6.13) u^^ u strongly in L''(51), as to — > 00, for any k : 1 < k < 7*. 
Since implies 2 < 7*, it follows from (pT^ that 

(6.14) — > u strongly in L^(il), as to — >■ cxd. 

Using the definition of (/. ()5.ip ) and the result (|6.14p . we can prove that 

(6.15) ® Ug^$e^ (|ue,^ I) — )■ u (g) u stroiigly in ($7) , as to oo. 

Then gathering the information of (|6.12p and (|6.15p . we see that G = u (g) u. 

Finally, using the convergence results (|6.9p - (|6.12p and observing (|6.15p . we can 
pass to the limit to — > oo in the following integral identity, which results from (j5.5p , 

^ (/.i|D(u,,jr-2 + /i2|D(u,„)rW-2) D(u,„) : D(^) dx 

- [ [ue„(gu,„,$,„(|ue„|) + F] :D((^)dx = 0, 
valid for all e V, to obtain 

(6.17) / (AiiSi+/i2S2-u®u-F) :D((^)dx = V S V. 

7. Determination of the pressure 

Since we shall use test functions which are not divergence free, we first have to 
determine the approximative pressure from the weak formulation (|6.16p . First, let 
w' be a fixed but arbitrary open bounded subset of such that 

(7.1) cj' CC fl and doj' is Lipschitz 

and let us set 
(7.2) 

:= (mi|D(u,„)P-2 ^ ^2|D(u,„)r('')-2) D(u,„J - ® u,„<i>,„ (|ue„ |) - F. 
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Using assumption ()4.3|) and the results (|6.4|) . (|6.6p . ()6.8p and ()6.10p . we can prove 
that 

(7.3) Qe,„ e L'-(r!), 
where r can be taken in such a way that 

(7.4) 1< r < ro := min |7', ^ 

Note that ro = min{7', 2(^3^} if 7 < ^ and we can take = 7' if > 7. Next, 
we define a linear functional 

(7.5) : Wyiu') ^ W-i''-(cj') 

by 

(7.6) (ne„,^)w-..(.oxw-'(.') ■■= I ^^r. ■■ D(^)dx. 

Using (|7.5p - (|7.6p . we can prove, owing to (|7.3p . that exists a positive constant C 
independent of m such that 

(7.7) ||n,,J|(v^,), <C7. 

Note that here V,.' is taken over w'. Moreover, since V is dense in V^', we can see, 
due to dnmi), dZll) and that 

(7.8) (n,„,¥.)(v^,),xv„, =0 V(^e V,-. 

By virtue of (|7.5p - (|7.8p and due to assumption (|7.ip . we can apply a version of de 
Rham's Theorem (/. Bogovskii [71 Theorems 1-4] and Pileckas [231 Section 1]) to 
prove the existence of a unique function 



(7.9) p,,„ e L'' (w'), with / pe,„rfx = 0, 
such that 

(7.10) (n,„,(p)^_,,,^^,^^^,,.,^^,^ = £p,„divvpdx V ^ e V^Y\u:') 
and 

(7.11) lbe,JlL.'(.0 < l|ne„||(V,,)'- 

Then, gathering the information of (|6.16p . (|7.2p . (|7.6p and (|7.10p . we obtain 



(7.12) 



f U|D(u,„)r-2 + ^2|D(u,„)rW-2\]3(^^^) . 

/ F:D((/3)dx+ / Ue„^ (g) Ue^$e„(|Ue„J) : D((/3) dx + / pe„^div(^dx 



for all e Wg'*^ (cj'). On the other hand, due to (|7.7p and (|7.1ip and by means of 
reflexivity, we get, passing to a subsequence, that 

(7.13) p^^ ^ po weakly in L'' (w'), as m — > 00. 
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Next, passing to the hmit m — > oo in the integral identity ()7.12p by using the 
convergence resuhs ()6.10p - ()6.12p . observing that in the last case, by virtue of (|6.15[1 . 
G = u (g) u, using also (|7.13p , we obtain 

(7.14) / (^iSi + ^^282 - u- F) : D(<^)dx = / padivipdx 

J u:' J u' 

for all ip e Wg'*^ (w'). Then, proceeding analogously as we did for ()7.5p - ()7.7p . we 
can define a linear functional 

(7.15) Ho : Wj''^'(w')^W-i''^(w') 

by 

(7.16) (Ho, <p)w-i.^(.')xwJ-'(..') - Qo ■ 

where Qo MiSi + ^282 — u (g) u — F. such that (|7.7p - (|7.8p are verified with Hq 
and Qo in the places of Ile^^ and Qe^ . In consequence, by the same version of de 
Rham's Theorem aforementioned, there exists a unique function 

(7.17) PoeL'''(w'), with f p^dx^O, 
such that 

(7.18) (no,^)^_,,(^,)^^,„.(^,^=£podiv^rfx, 

and (|7.1ip is verified with and Ho in the places of pe„, and Ile^ . Then gathering 
(fTTB)) and (TTTSI . we achieve to 

(7.19) / (/^iSi + ^282 - u® u - F) : D(.^) dx = / p^divifidx 

J uj' J u' 

for all (f S Wq''^ (w'). Finally, combining (|7.14p and (|7.19p . and by means of 
uniqueness, we conclude that 

8. Decomposition of the pressure. 

The main idea in this section is the application of a method to locally decompose 
the pressure found in the previous section. For that, we shall use a lemma which 
is proved by using a direct decomposition of . which in turn is equivalent to the 
weak L^-solvability of the Dirichlet problem for the Bilaplacian in bounded domains 
with boundaries (/. Simader and Sohr [SS]). With this in mind, let w be a fixed 
but arbitrary domain such that 

(8.1) ojCCuj'ccfl and 9a; is C^. 

Lemma 8.1. Let 1 < s < 00. Then for every v* G W^^'^(a;) there exists a unique 
u G WQ''^(aj) such that 

AuA^bdx = (t'*,0)w-,=Mxw„^-'(..) ^ <^ e C^iu;). 
Proof. See e.g. Wolf [SI Lemma 2.4]. □ 
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To simplify the notation in the sequel, let us set 

A'ioj) := {a e L^(w) : a = Am, ue W^' "(w)}. 
Applying Lemma 18. 1[ with 5 = 7' first and then with s 



1 < s < 00. 

attending to the 



definitions of A''' (w) and A 2 (w), and using (|6.10p and (|6.f f p by one hand and 
(|6.12p and (|6.f 5p on the other, we can infer that exist unique functions 

(8.2) pi^GA^'H, 

(8.3) p^^eA^H 

such that 
(8.4) 

J UJ 

/ii (|D(u,„)r-2D(u,„) - Si) + /i2 (|D(ue„)r(^)-2DK„) - S2 



(8.5) 



/ pl^^A(j>dx = (U£^ (g) Ue 

J UJ J UJ 



u (X) u) : V^0cJx 



for all (f) e Cg°(aj). The next result helps us with the estimates of the local pressures 
found in (lOl-dS^l). 



Lemma 8.2. Let 1 < s < 00 and assume that p G L''(ti;) and Q G L*(aj) are such 
that 



f pA(t)dx= f Q:V^(j)dx V0eCS°(a;) 

J UJ J UJ 



The 



IWIl^M <C||Q||l.m, 
where C is a positive constant depending on s, N , uj and on the Calderdn-Zigmund 
inequality's constant. 

Proof. See Wolf [SlJ Lemma 2.3]. □ 

Attending to (|6TT|) . ((6T0| and jS^) by one hand, and ((6T2| . ((6T5)) and jHS]) on 
the other, a direct application of Lemma 18.21 with 5 = 7' and then with s = 
yields 

IbLllL.'M <Ci|||D(u,„)r-2D(u,„) - Si||Ly(^) + 
C2|||D(u,„)rW-2D(u„J-S2||Ly(^), 



(8.6) 



(8.7) 



I XI, < C3||Ue„^ (g)Ue„$£„(|Ue„|) - U®U|| ^ 
L 2 (lj) L 2 [uj) 



where Ci, C2 and C3 are positive constants depending on 7, iV, w and on the 
Calderon-Zigmund inequality's constant. 

On the other hand, combining ()7.12p and (|7.14p . and using the definition of the 
distributive derivative, we obtain 



div (|D(u,„)r-2D(u,„J - Si + |D(u,,jrW-2D(u,„J - S, 
div(u,„ (gu,„$,^(|u,„|) -u(8)u) = V(p,„ - po) 



in V'{uj). 
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Then, testing ()8.8p by V<p, with e C3°(cj), integrating over u and comparing the 
resulting equation with the one resulting from adding (|8.4p and ()8.5|1 . we obtain 



Inserting this into (|8.8p . it follows that 
(8.9) 

div (|D(u,,jr-2j3(^j^^) _ ^ |D(u,„jr('')-2D(u,„) - S2) - 
div(u,„ ® U£^$e„(|u<:„J) - u(8)u) = V +pL) 



in V'{uj). 



9. The Lipschitz truncation 
To start this section, let us set 
(9-1) := (u£„^ - u)xtj, 

where denotes the characteristic function of the set w introduced in (j8.ip . Having 
in mind the extension of (|8.9p to R^, here we shall consider that 

(9.2) T,„^ T,i^^ + 

is extended from w to by zero, where 
(9.3) 

Ti := - f|D(ue„)r-2D(u,„J - Si + |D(u,,jrW-2j3(^^^j _ g A ^^1 



(9-4) Tf^^ := Ue„ «) U<,^$e„(|Ue„|) - U(g) U + pf^J, 

and I denotes the identity tensor. 

Now, due to the definition (|9.ip and by virtue of (|6.9p and (|6.13p . we have 

(9.5) w,,„^0 weakly in Wi''^(R^), asm ^00, 

(9.6) Wg,^ — > strongly in L''(R^), as m — > 00, for any k : 1 < k < 7*. 

Moreover, due to ((57TU)) . and ([SIB)) by one hand, and due to ((O^ . 

and (|8.7p on the other, we have 

(9.7) l|TfLllLV(E«) 

(9.8) II ^ <C. 

In addition to (j9.8p . we see that, due to (|6.13p and (|8.7p . 

(9.9) T2^->0 strongly in Lt(R^), as m — > 00, for any k : 1 < < 7*- 

Next, let us consider the Hardy-Littlewood maximal functions of Iw^^^ | and 
|Vwe^| defined by 

>'(|we„|)(x) sup I [ |w,„(y)|dy, 

X(|Vw,,J)(x) := sup — — 1-— /■ |Vw,„(y)|dy; 

0<i?<oo i-N[-DB(X)) Jbr{x) 

where i?i?(x) denotes the ball of R^ centered at x and with radius R > 0, and 
£Ar(a;) is the iV-dimensional Lebesgue measure of w. Arguing as in Diening et 
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al. \12\ p. 218] and using the boundedness of the Hardy-Littlewood maximal oper- 
ator Ai (see e.g. Stein |30l Theorem 1. 1.1]), we can prove that for aU 7ti G IN and 
all j £ ¥(q there exists 

(9.10) X,^,g[2^\2''^') 
such that 

(9.11) £A,(i^™,,) < 2-^A-^|lwe„|lL„(R«), foranyK:l<At<7*, 

(9.12) (G„,,) < ||Vw,„||l.(i,«), 
where 

F,„,, := {x e : A^(|w,,J)(x) > 2A„j} , 

Gra,j {x e : X(|Vw,^|)(x) > 2A^,,} . 

Setting 

(9.13) Rm,j '■= U Gm.j U {x G : X is not a Lebesgue point of |we^|} , 
we can see that, by virtue of ^JT^i-^J^ and 



(9.14) limsup/Iiv (i?mj) < limsupC2"^A 



-7 



The following result helps us to approximate M^^''''-functions by Lipschitz ones. 

Lemma 9.1. Let cu C M,^ be an open bounded set with a Lipschitz- continuous 
boundary dui. Assume that s > 1 and let u G W^'*(cj). Then, for every Ai, A2 > 0, 
there exists v G W"'"'°°(w) such that 

I|v|1l-(c^) < Al, 

I1Vv|1l=o(„) < CA2, C^C{N,Lo). 

In addition, 

{x e w : v(x) ^ u(x)} C 

n ({x e M,^ : 7W(|u|)(x) > Ai} U {x e : Al(|Vu|)(x) > AajU 
{x G R^ : X is not a Lebesgue point of |u|}) . 

Proof. See Acerbi and Fusco |T] for the original result and Landes |181 Proposition 
2.2] for the last statement. □ 

Then, by Lemma [01 together with the definition of w^^ (/. (|9.ip ). there exists 

(9.15) z,„j e W ' (R ), z„j = I Q R^\a; ' 
where 

(9.16) := {x e w : z,„j(x) ^ We„^(x)}, 
such that 

(9.17) |lZmj|jL=c(^) < 2Amj, 

(9.18) ||Vz™,,|1l=o(^) <CA™,„ C = C{N,oj). 
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Moreover, by the last statement of the above lemma and using the notations (|9.1ip - 
([9l3| and ([QAB]) . 

(9.19) A,„„j di^nRrnJ. 

As a consequence of (|9.14p and (j9.19p . 

(9.20) lim sup Cn {A^j) < lim sup C2-^ . 

m— >-oo m— >-oo 

On the other hand, due to (PT7| - (P7T5)l and we can prove that for any 

J e Mo 

(9.21) Zmj^O weakly in Wj'''(a;), as m — )- oo. 
Then by Sobolev's compact imbedding theorem, we get for any j £ INq 

'^m.j strongly in 'L'^{lu), as m — > oo, for any k : 1 < k < 7*. 
Using this information, (|9.17p and interpolation, we prove that for any j G INq 

(9.22) Znij strongly in L'*(a;), as m 00, for any s : 1 < s < co. 
Finally, as a consequence of ()9.2ip and (|9.22p , we obtain for any j € 

(9.23) z,„j weakly in Wo''*(cj), as m — > 00, for any s : 1 < s < 00. 

10. Convergence of the approximated extra stress tensor 

Let us first observe that, using the notations (|9.2p - (|9.4p . we can write (|8.9p in 
the following simplified form 

(10.1) divY,„ = in V'{uj). 

On the other hand, due to (|9l2l)-(|9ll), Y,„. e L''(K.^) for r satisfying to (fr4| . 
Then, using this information and ()9.23p . we infer, from ()10.ip . that for any j £ INq 

(10.2) I : Vz,„j dx = 0. 

Expanding T^^ in p0.2p through the notations (|9.2p - (|9.4p and subtracting and 
adding the integral 

f (Mi|D(u)r-2 + Ai2|D(u)r(-)-2) D(u) : D(z„,,) dx 
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to the left hand side of the resulting equation, we obtain for any j e 

ti2 [ (|D(u,„)|''(^)-2d(u,„)-|D(u)|«(^)-2d(u)) :D(z„,,)dx = 

J UJ 

//I f (Si-|D(u)r2D(u)) :D(z™.,)dx+ 
(10-3) f,, f (S2 - |D(u)|«('')-2D(u)) : D(z™,,)dx+ 

J UJ 

J CO 

/ (u,„ ® Ue„<I>e,„(|ue„|) - u® u + p^^l) :D(z™j)dx 

- /I +72 . + +7'^ . 
We claim that, for a fixed j, 
(10.4) 

/.I lim /" (|D(u,„jr-2]3(^^^j_|j3(^)p-2]3(^)^ . 

m— >oo / 

lim / (|D(u,„)|''(^)-2d(u,„J-|D(u)|'?W-2j3( A.j3^^ .^^^<^2-i 

To prove this, we will carry out the passage to the limit m — > 00 in all absolute 
values I jl, i = 1, . . . ,4. 

• limsup„_^o2 \ Jl^ j l ==0. By (|9.23p . with s = 7, this is true once we can justify 
that Si — |D(u)|'''^^D(u) is uniformly bounded in L'*' (w). But this is an immediate 
consequence of (|6.10p . 

• limsup^_^oQ \ J'^n j\ — ^- Again, by (|9.23p with 5 = 7, this is true if S2 — 
|D(u)|^^''-*~^D(u) is uniformly bounded in (w). Analogously to the previous 
case, we can justify this now by using (|6.1ip . 

• limsup^^oQ \ Jm.j\ ^ C2~^ . In fact, by Holder's inequality and (|8.6p together 
with (|6.10p and (|6.1ip . and using the definition of (/• (|9.15p ) together with 
the fact that divw^^ = in w (/. (P7T|) ). 

< C'i||divZ™j|jL,3(<^) < Ci|lVz™j|lL-(c^)£Ar(A,„j)T. 

Then, by the application of ((^1^ and (|O0)) . it follows 

(10.5) \Jl,.^ < CalimsupA™,, {2'^ X'^)- < ^22"^. 

m— f 00 

• limsup,^_^Q^ \Jmj\ — ^- Using Holder's inequality and the notation (|9.4p . we 
have 

\JmJ < l|T"f„J|Li(t^)||Vz,„.j||L==(c^) 

< CillT^^JlLi(t^) 0, asm^cx). 

The last inequality and the conclusion follow, respectively, from (|9.18p and (|9.9p 
with K = 2, observing that here the assumption (|4.2p implies 2 < 7*. 
Gathering the estimates above we just have proven (|10.4p . 
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We proceed with the proof by using an argument due to Dal Maso and Murat |9l 
Theorem 5] (see also the references cited in [9l Remark 4]). Firstly, observing the 
definition of Zmj (/• (|9.15p ). we have 

Ml f {\T>{u,J\'^-^-D{u,J~\-D{u)p~^-D{u)) ■.■D{z^,,)d^+ 

J UJ 

(If^-^) M2 / (|D(u,„)r(-)-2D(u,„) - |D(u)r(-)-2D(u)) : D(z,„,,) dx 

J u 

where 

II^ := / (|D(u,„jr-2D(u,„J - \T>iu)r'T>iu)) : (D(u,„) - D(u))dx, 
Il.j (|D(u,„jr(-)~2D(u„J - |D(u)r(-)-2D(u)) : (DK„) - D(u)) dx, 

III^^ / (|D(u,„)r-2D(ue„) - |D(u)r-2D(u)) : D(z,„,,)dx, 
III^ := 1^ (|D(u,„)rW-2D(u„J - |D(u)rW-2D(u)) : D(z,„,,)dx. 
Then ((To:4)) and ([TOJl imply that 

(f 0.7) lim sup (/i,,^. + /2„ < C lim sup (| J/,^^ + /I,^,, I) + 2"^ 

For the terms 11^ ^ and //^ j , we have by applying successively Holder's inequality, 
([6l0| and ((6TT|l . altogether with ([9T^ . 

Cili (|D(u,,jr-2D(u,,J - |D(u)r-2D(u)) |Il.'(a,„.,)I|Vz,„,,||l.(a„,,) + 

(|D(u,„jr(-)-2D(u,„) - |D(u)rW-2D(u)) |lLy(A„„)l|Vz™,,||L.(A„„) < 

Then using the same arguing as we did for (jlO.Sp it holds for any j S 

(10.8) lim sup + \ < . 

As a consequence of (|10.7p and pO.Sp , we obtain for any j e INq 

(10.9) limsup + < . 

Now, by Holder's inequality and having in mind the definition of z^j (/• (jQ.lSp l. 
we have for any 9 G (0,1) 

g!^ dx <Ci [iliJCNito \ A„.,)i-^ + (IT,i„^^.)'^iv(Anj)'"'] + 
[{lUCNiu;\A^,,y~'' + (II^^^^O'^C^lAn,,)'"'] 



(10.10) ^ 
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where 



.Ml ||D(u,„)r-2D(u,„) - |D(u)r-2D(u) 



M2 



and 11]^ j, 1-m.j ^^^d are, respectively, the functions I^^^ j, ^^m,ji ^mj 
and J, but with the integrand functions replaced by their absolute values. Ar- 
guing as we did to prove p^?5)) - (fTU?g)l and using ([11201), it follows from pO.lOp 
that 



(10.11) limsup / ,gf^^ dx < Ci2~^t + C22' 

m— >-00 J CO 



Since 7 > 1, G (0, 1) and j e INq is arbitrary, 2^"^ and 2-^t-(i-^)J ^ 0, as 
j 00. This and pO.lip imply that for any 9 £ (0, 1) 



limsup / g^^ dx = 0. 

m— ^00 J uj 

Then, passing to a subsequence, 

(10.12) g^^ — > a.e. in u, as m — s> 00. 

Due to the fact that |D(u)|'''~^D(u) and |D(u)|*('')~^D(u) are strictly monotonous 
and continuous tensors on D(u), we can apply O Lemma 6] together with (|10.12p 
(see also Lions [191 Lemma 2.2.2]), to establish that 

(10.13) D(ue,^ ) D(u) a.e. muj, as m — > 00. 

Finally, ((0)l - (p3)l and (|10.13p allow us to use Vitah's theorem together with ((OU)) - 
mJT\i to conclude that Si = |D(u)|''-2D(u) and S2 = |D(u)|'?(^)-2j3(„)^ 
This concludes the proof of Theorem 14.21 ■ 



Remark 10.1. Using the parabolic versions of the techniques we have used here 
(see e.g. Oliveira |21| ). we can establish an existence result for the parabolic model 
of the problem lll.S^) - ![T75\) . This work will be the aim of a forthcoming paper. 
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